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I. INTRODUCTION 

The studies of the explosive processes like SNe and 
GRBs have been revealing the evidences of compact ob- 
jects such as neutron star and/or black hole at the cen- 
ter of explosive events. And it is likely that massive 
black holes are responsible for the activities of AGN. 
The existence of a black hole leads to the idea that the 
energy sources for powering violent astrophysical phe- 
nomena like AGN and GRB are likely to be the grav- 
itational binding energy of the accreting material onto 
a black hole and the rotational energy of a black hole 
itself P, 1^ It has been also expected that the 

strong and large scale magnetic field configuration can be 
formed around the compact objects at the centers, par- 
ticularly for the late stage of the core-collapses of rapidly 
rotating and strongly magnetized stars. It is then nat- 
ural to suppose a simple physical situation in which the 
energy of the accreting material and the rotational en- 
ergy of a black hole can be extracted out through the 
magnetic field lines in the form of the Poynting flux. 

The Poynting flux has been considered to be a viable 
route to transport the energy and angular momentum 
along the magnetic field lines and one of the efficient ways 
to tap the black hole's rotational energy using magnetic 
field lines anchored on the horizon 1] supported by the 
external current. The relevance of using the Poynting 
fiux in describing the powerful and highly collimated as- 
trophysical jets observed in AGN and quasars has been 
suggested long time ago ^'^d many interesting works 
have been developed. One of the characteristics of the 
Poynting flux is that it carries very little baryonic com- 
ponent compared to the hydrodynamic flow. This prop- 
erty of Poynting flux is found to be consistent with the 
required property for powering GRB Q . Recently numer- 
ical simulations also have begun to demonstrate the elec- 
tromagnetic energy extraction 6, X 8] from a black hole 
surrounded by the force-free magnetosphere, which can 
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be established when the electromagnetic field is strong 
enough with sufficient charged particles for space-charges 
and currents and the inertia of the plasma can be ig- 
nored. The Poynting flux in a system of black hole- 
accretion disk has also been studied in connection with 
GRBsHlllIll and it is found that the evolution of the 
system largely depends on the Poynting outflow from the 
diskjlllllllllll. 

In this work, the nature of the Poynting flux from the 
horizon and the ergosphere of a black hole and from the 
accreting flow onto a black hole will be discussed. To 
make the discussion simpler the environment through 
which the Poynting flux carries energy out is assumed 
to be force-free with steady and axisymmetric conflgura- 
tion and the accreting flow is assumed to be conflned on 
the equatorial plane. In section II, the force-free magne- 
tosphere is introduced together with the consistent con- 
straint and the basic features of the Poynting flux is dis- 
cussed in section III. In section IV, the accretion flow 
dominated by the Poynting flux is discussed briefly and 
the possible Poynting flux from the equatorial plane in- 
side the ergosphere is discussed in section V. 

II. FORCE-FREE MAGNETOSPHERE IN KERR 
BACKGROUND 

The force-free magnetosphere can be established when 
the electromagnetic field is strong enough with sufficient 
charged particles for space-charges and currents and the 
inertia of the plasma can be ignored. The role of the 
plasma in the force-free magnetosphere is to provide the 
charge and current sources for the field. The mathemat- 
ical expression for the force-free limit is given by 

Ff,,r^o, (1) 

where J" is the electromagnetic current density. In terms 
of the physical quantity defined by FIDO, it can be writ- 
ten as 

pE + JxB^O, (2) 
from which one can also observe that 

E-B^O, (3) 
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which imphes that the electric potential is conserved 
along the magnetic field lines. It should be noted that 
the force-free condition is not identical to the ideal MHD 
condition, 



or equivalently 



F^.u" = 0, 



E + VX B ^0, 



(4) 



(5) 



where u'^ is the four velocity of the plasma. 

The force-free requirement should be consistent with 
the condition that there is no local frame in which a pure 
electric field can be seen[l6l \v^ . In other words, the 
invariant 



- 



^F^'^F^, >0, 



(6) 



has to be always checked. On the other hand, when one 
can find a region where Eq.® can not be consistent it 
is very likely that it is the signature of a break-down of 
force- free condition. Hence it can be a very useful method 
to find out the non-force-free region (most likely a part of 
the region) even without knowing the detailed knowledge 
of the plasma around. 

In the steady and axisymmetric case around the Kerr 
black hole, we get from the the force-free condition, 
Eq.O, 



E = [VLf +/3) e"^ X BP. 

a 



(7) 



f2i?, which is defined by dA^ — —U,pdA^^ is the angular 
velocity of the magnetic surface that rotates rigidly in an 
axisymetric and stationary state Throughout this 

work, the background metric is assumed to be the Kerr 
metric 10. Using the Boyer-Lindquist coordinates 20J in 
the natural unit G = c = 1, the non vanishing compo- 
nents are given by 



500 



-{a -w (3 ), = 500 = /3 

9r 



A 



7 gee = P^, 94,4, = w^, (8) 



where 



a 



A = 



pVA 



2aMr E 
- 2M r, = + cos^ 
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2j = (r +a ) — a Asm 



(9) 



Using Eq. ITJ , the invariant B^ — E"^ can be written as 

B^-E^ = -f{nF,r,0){BP'f + {B^f, (10) 

where B^ is the poloidal component of the magnetic field 
and 

f[VLF,r,e) = ^\gci>ci> + '2.^Fgo4> + 9oa 



(11) 



As a simple example, consider a region very near to the 
horizon, a ^ 0. Then Ea. HlU|) can be written as 



B^ - E^ 



H) 



^2^ 



{BPf + (5-^)2.(12) 



One can see that without the toroidal component of the 
magnetic field, B'^ , the above equation becomes negative 
and the force-free condition breaks down. 

The force-free condition is assumed for the plasma with 
a strong magnetic field, which is sufficiently tenuous that 
the plasma exerts no force on the magnetic field. If the 
whole region is force-free, there is no place to extract 
energy for the Poynting flux. We suppose a Poynting 
flux as a viable route to transport the energy and an- 
gular momentum along the magnetic field lines. Then 
there should be a region where matters exerts forces on 
the magnetic field such that they convey the energy and 
angular momentum along the field lines in the form of 
Poynting flux. For the Poynting flux out of accretion 
disk, it is naturally expected that the force-free condi- 
tion holds only for the magnetosphere not on the disk. 



III. BASIC FEATURES OF POYNTING FLUX 

To discuss the basic feature of the Poynting flux, sup- 
pose an artificial surface at large radius which encloses 
the black hole and accretion disk for example. The elec- 
tromagnetic field near the surface is described by the the 
Maxwell equation, 



' ■ 7/ 



47rJ^. 



(13) 



The conservation of bulk current density should be mod- 
ified since the bulk current density has a discontinuity on 
this arbitrary surface|21j. To ensure the conservation of 
the current with this boundary, the total conserved cur- 
rent density should be the sum of the bulk current 
density and the surface current density j^: 



(14) 



where the effective surface current density are determined 
using the current conservation law 



0. 



(15) 



As an example, let us consider a spherical surface for 
an axisymmetric and steady case. One can have 

If the radius of the arbitrary surface is located at far 
enough distance from the central object, then we can 
use a flat space-time geometry. Then the energy flux of 
electromagnetic field is given by 



1 

47r' 



.E^B'f 



(17) 
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which is the radial component of the Poynting flux. On 
the arbitrary surface using Eq. H16|) , the Poynting flux can 
be written as 



(18) 



One can note that for a non- vanishing Poynting flux the 
force-free condition cannot be maintained on the bound- 
ary, j ■ E ^ 0. It corresponds to the 'effective battery if 
it is positive (or 'ohmic dissipator' if negative). Similarly 
the z-component of the angular momentum flux can be 
written by 



(19) 



where Ea. ((TC|l is used and zu is a. cylindrical radius. It 
corresponds to the effective 'torque' in z-direction ex- 
erted by the surface current although it is force-free in 
the original setting. This observation implies that the 
non-force-free nature of the system sitting at the cen- 
ter inside the surface can be manifested on the arbitrary 
boundary surface outside. 

Now put the arbitrary surface near to the gravitating 
object with the smallest radius but big enough to enclose 
the energy source while keeping outside to be force-free: 



F^,r = 0. 



(20) 



We need a fully relativistic treatment. For Kerr geome- 
try, the energy fluxj2^ is given by 



(21) 



Using Eq.(|Zj), for a rigidly rotating magnetic field line 
with angular velocity Qp in a, force-free environment up 
to the boundary surface, Ea. H21|) can be written by 



47r (f^F + /3) 



(22) 



IV. POYNTING FLUX DOMINATED 
ACCRETION FLOW 



We consider a system of a rotating black hole and a 
two-dimensional accretion disk, which is assumed to be 
on the equatorial plane in the background metric of a 
Kerr black hole and the accretion is supposed to be driven 
by the Poynting flux,23j 

The stress-energy tensor is decomposed into two parts, 
the matter part, Tj^'^ , and the electromagnetic part. 



(24) 



In general the stress-energy tensor for the accretion disk 
is determined by mass density, internal energy, pressure, 
viscosity, radiative transfer and etc. (2^. Since we are 
interested in the accretion dominated by the Poynting 
flux, the two-dimensional accretion disk is assumed to be 
non- viscous, cool, and non-radiative such that the matter 
part is given by 



(25) 



where pm is the rest-mass density and w'' is the four ve- 
locity of the accreting matter. It is also assumed that 
there is a negligible mass flow in the direction perpendic- 
ular to the disk: — 0. We can obtain the energy flux 
and the angular momentum flux given by 



CM 



(26) 



For an idealized thin disk on the two-dimensional plane 
we take 



p^^^S{9-n/2), 
P 



(27) 



where am is the surface rest-mass density. The rate of 
the rest-mass flow crossing the circle of radius r defines 
the mass accretion rate M+ by 



which shows a typical form of Poynting flux. 

Using the effective currents defined above we can ob- 
tain 



-2TramPU^- 



(28) 



(23) 



We can note that the non-force-free nature of the source 
inside for the Poynting flux outside is represented by non- 
vanishing which 'flows' across the magnetic field lines 
on the boundary surface. Basically if the whole region 
is force-free then there is no way to get a Poynting flux: 
no energy source. Hence the magnetic field lines should 
be anchored onto a non-force-free region from which the 
Poynting flux can carry out the energy. However the 
detailed nature of non-force-free region can only be un- 
derstood in the frame work of the relativistic magneto- 
hydrodynamics. 



It is identical to the one derived in j2J| and [23 , in which 
the vertical structures of the disk are integrated. The sta- 
tionary accretion flow implies that is r- independent. 
Using Ea. H26|) the radial flow of the energy and an- 
gular momentum of the matter at r can be given by 
MoM+ and — Af+ respectively. 

The horizon of a black hole is a mathematically well 
defined surface, on which the appropriate surface current 
density has been studied in depthUyj . On the other hand, 
the surface of the physical accretion disk is not expected 
to have any sharp boundary and various shapes have 
been suggested depending on the properties of the accre- 
tion flow. In this work, however, we assume a simplified 
accretion disk, which is vertically squeezed down to the 
equatorial plane, a two-dimensional accretion disk. Then 
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we can associate the surface current density in Ea. H14() 
in a simple way, 



1 



(29) 



which are nothing but Gauss' law and Ampere's law as 
given by 



1 



-1-Bi K^^l-Bl (30) 



where (Tg and are surface charge and surface current 
density(spatial) respectively and the reflection symme- 
try with respect to the equatorial plane, F^'^ = —F^''\ 
is used. The surface current densities are responsible for 
the Poynting flux from the disk. Hence one can consider 
the singular equatorial plane with non- vanishing surface 
currents manifests the non force-free nature of the accre- 
tion flow driven by Poynting flux. 

The energy and the angular momentum of the disk can 
be carried out by Poynting flux along the magnetic field 
lines which are anchored on the disk. In our simplified 
model, it is the main driving force for the accretion flow. 

Using the Killing vector in i-direction, ^'^ = 
(1, 0, 0, 0), we can define the energy flux from the 
stress-energy tensor of the electromagnetic field, T^'^, 



'-'EM 



where 



47r 



(31) 



(32) 



For the steady and axisymmetric case in this work, 
E"^ = 0, and we get the total energy flux of electromag- 
netic field in ^-direction, Sf^, given by 



D 



P 



(33) 



The first term in Ea. (|33|l can be considered to be an 
ohmic interaction, E ■ K. For the current in the direction 
of the electric field tangential to the disk there might be 
an energy dissipation into the disk surface. This is what 
one can expect on the black hole horizon J^]. However 
for the current in the opposite direction this term cor- 
responds to the electro-motive force and it is the case 
for the accretion disk on the equatorial plane discussed 
by Blandford and Znajek0, Qas well as in this work. 
The second term can be interpreted as a magnetic brak- 
ing power on the rotating body with an angular velocity 
LJ = -13. 

Using the Killing vector in (/)-direction for the axial 
symmetric case, 77^ = (0, 0, 0, 1), we can get also the 
total flux of the angular momentum in ^-direction, 



C 



D 



2ttp + + 



1 



wK'Bl, 



(34) 



which is nothing but a magnetic torque exerted on the 
surface current density K^. 

Now the Poynting power measured at inflnity is given 
by 

Pd = - J drr{-aElB'l+ I3wbIbI), (35) 

and the rate of angular momentum transfer, Pj^, out of 
the accretion flow is given by 



drrwBlB\. 



(36) 



A. Accretion Equations and Stream Equation 

The dynamics of an accretion flow is determined by 
the conservation equation of the stress-energy tensor: 



(37) 



From the conservation of stress-energy tensor, one can 
obtain the accretion equations 126*1 driven by the Poynt- 
ing flux in a force-free magnetoshpere given by 



{drUo)M+ + rmVLpBlBl = 0, 
{drU^)M+ - rwBlB\_ = 0, 

drUQ 



g —[drU^) 



27rr2 
/A 



2'Kr' 



BiBi \ —nl - 1 ) = 0, 



(38) 
(39) 



(40) 



where the angular velocity of the disk is given hy fin = 
u'^/u'^. Eas. (|38|) and H39|l correspond to the energy and 
angular momentum conservation respectively in the sta- 
tionary and axisymmetric setting in this work. In Ea. l|38|) 
the radial change of the energy flux of accreting mat- 
ter, —U0M4. is balanced by the Poynting flux in theta 
direction [231, Sf)- Similarly the radial change of the an- 
gular momentum flux, u^M^ is balanced by >C^. The 
radial equation Eq. (|40|) essentially determines the orbital 
motion of the disk. In the absence of external fields in 
Ea. (|40() . the angular velocity is determined as 



drUo 



(41) 



which is one of the characteristics of the Keplerian orbit. 
Hence deviations of from the Keplerian one are natu- 
rally expected in the Poynting flux dominated accretion 
disk. 

It should be noted the force-free condition holds only 
for the magnetosphere not on the disk. In fact the driving 
force for an accretion flow on the disk is the magnetic 
braking in Eqs.(|2HI) andllSHl- Then Eq.(|2nil becomes 



Pjj = / drr(voQ.FBlBl). 



(42) 
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From Eas. (|38|l and (|39() . one obtains an interesting re- 
lation 



drUQ 
drU^ 



(43) 



which has no explicit dependence on the field configura- 
tion. It seems to imply that ilj? is determined essentially 
by the dynamics of the disk. However one should note 
that the dynamics itself is governed not only by the grav- 
ity but also by the electromagnetic field as well. 

The configuration of the ordered magnetic field around 
the accretion disk and the Poynting outflow from the 
disk has been discussed both in analytical and numerical 
studiesj27t .28. 29. .30] . For the strong enough electro- 
magnetic field around the compact object, the force-free 
magnetosphere can be established. In the non-relativistic 
formulation, Blandfordjj] suggested an axisymmetric and 
stationary electromagnetic field configuration around an 
accretion disk on the equatorial plane. The poloidal field 
configuration for a black hole is known to satisfy a second 
order elliptic differential equation called Grad-Shafranov 
eauation|i31l| or a stream eauation[32i| for the stream func- 
tion ^ and current /. The poloidal and toroidal compo- 
nents of the magnetic field can be written in terms of ^' 
and / respectively: 



1 



2nm 



21 

tua 



(44) 



Possible types of solutions in a forcc-frcc magneto- 
sphere has been discussed recently in the relativistic 
formulation [s^. Is^. To calculate the Poynting flux we 
need to know the conflguration of the magnetic held 
which is consistent with the accretion equations. And 
the problem is reduced to solve the four coupled equa- 
tions, three accretion equations and a Grad-Shafranov 
equation(stream equation). 



B. Example: Numerical Solutions in Schwarzschild 
Background 



of the solutions, ^q, suggested by Blandford and Znajek 
is given by 



^-0 = nCX, 



(46) 



where 



X = r(l =F cos 61) + 2Af (1 ± cos6l){l - log(l ± cos6l)X.47) 

For non- vanishing flp and /, we consider a solution for 
which the shape of the magnetic surface are the same as 
the magnetic surface defined by Vl/oH^- Thus 5',/ and 
rip are assumed to depend on X. We suppose that the 
derivative of has the same form as in the flat back- 
ground, Ea. (|46|l . That is, we take an Ansatz such that 



dX 



ttC 



(1 



r!2,X2)l/2- 



(48) 



Four basic equations governing the accretion flow un- 
der the infiuence of paraboloidal-type configuration are 
solved numerically [35j. The radial variations of uq, u^, 
I and Up on the disk are basically functions of the ac- 
cretion rate and the strength of the magnetic field. Nu- 
merical calculations show that the angular velocities of 
the magnetic field lines (Oi?) are different from either 
the Keplerian angular velocity (fi^) or the disk angular 
velocity(ri£)). It implies that the two-dimensional ap- 
proximation of the perfect conducting disk, for which 
rip ~ Qp), rnay not be implemented particularly with 
the paraboloidal type configuration. As expected in the 
previous section for the velocity of the magnetic field line 
vp{= rflp) less than the speed of light. Hp is found to 
be larger than Up). 

It is also found that the strength of the magnetic field 
is increasing as r goes near the inner edge. Since the mag- 
netic field as well as the surface current are increasing as 
r gets smaller, it is naturally expected that the Poynting 
fiux increases substantially as r approaches to the center 
. Although the numerical calculation does not go beyond 
r < 6M, it may indicate a possible electromagnetic jet 
structure near the inner edge of the disk. 



The coupled differential equations in a Kerr geometry 
in the previous section is not easy to solve. To show some 
examples of the solutions we try a numerical solution '35] 
in a simpler background, Schwarzschild background. In 
the Schwarzschild background, the stream equation is 
given by 



T V sm u 



f^F sin^ Qdr {r^VLpdr^) - 



V. POYNTING FLUX FROM ERGOSPHERE 

On the horizon inside the ergosphere, we can make use 
of the Znajck's boundary conditionfl^. 



In general there is no known analytic forms for the solu- 
tions. In the limiting case when rip and / vanishes, one 



and we get 



amOde {smeflpde'i') 
(45) 



- 



{By 



(49) 



(50) 



It implies that for the magnetosphere with = on the 
horizon the force-free nature can be justified and the re- 
gion where the electromagnetic field lines can exert non- 
zero force should be inside horizon. As discussed in sec- 
tion III, we can define a non-zero surface current on the 
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horizon boundary such that the rotation of black hole can 
be slowed down or speeded up and we get the Poynting 
fluxP|. Using Eg. 1)22(1 . the rate of energy extraction at 
infinity can be evaluated using the electromagnetic field 
on the horizon, 

Phole = ~ [ ^^F(SlF-^HW{Bydedcj). (51) 
Jr=rH 

It is to be noted that the force-free condition might not 
be maintained on the horizon for ^ 0. It may corre- 
sponds to the case when the force- free condition is invalid 
in the vicinity of horizon ,39j. which, however, depends on 
the details of the physical properties of the plasma. 

We will explore whether similar analysis can be applied 
to the magnetic field lines threading on the equatorial 
plane inside the ergosphere^^. The overall environment 
and physical processes are assumed to be axisymmetric 
and steady in the background of Kerr geometry. The ef- 
fects on the geometry due to the energy density of the 
electromagnetic field and the accreting material are as- 
sumed to be negligible. In this work we assume a very 
thin accretion flow on the equatorial plane. Let us sup- 
pose the force- free relation, F^" Jj^ = 0, holds up to the 
equatorial plane to see under what condition — E^ > 
is violated. Because of the reflection symmetry in the 
equatorial plane, there is no current crossing the equato- 
rial plane and B'^ = on the equatorial plane. Then we 
get 

B^-E' = -f{nF,r,e = 7T/2)^-^, (52) 

where f — determines the light surfaces [ioj. For a 
given r, the sign of B^ — E^ is then determined by /. 
Hence the condition for force- free requirement , / < 0, 
determines the acceptable range of the angular velocity, 
rii?(f), of the field lines which cross the equatorial plane 
at r. One obtains 

n-{r) < Qpir) < n+ir), (53) 

where 

n± = QzAMO ± — , ^ZAMO = (54) 

For such magnetic field lines which rotates with (r) < 
57(r) < f2_|_(r), we cannot expect any Poynting flux since 
it is simply force free. In other words, we cannot get 
non- vanishing B"^ for the Poynting flux. It is because for 
the force-free magnetosphere, Bt = waB'^ is constant 
on the magnetic surface and therefore B"^ should remain 
to be zero not only on the equatorial plane but also off 
equatorial plane along the field line. 

However for the field line with the angular velocity 
smaller than ri_, the region on the equatorial plane can- 
not be consistent with the force- free requirement. Hence 
the magnetic field line with U,p < VL_ can exert a force on 
the non-force-free region near the plane. And we can get 



a Poynting flux along such fleld lines, which pass through 
the non-force-free region. Consider a non force-free re- 
gion with finite width around equatorial plane. Since it 
is not a force-free region, the toroidal component can be 
developed along the field line off the plane, although on 
the equator it should be zero. It reaches the boundary of 
force- free region with B*^ = _B^, which is determined by 

-f{^,r,e)^-^ + {Bif^Q. (55) 

We can guess is not much different form fi_: 

nF = n^+5, (56) 

where 5 is determined by B"^. It also implies that the 
Poynting flux out of equator is possible only inside the 
ergosphere since f2_ = on static limit and outside the 
ergosphere there is no Poynting flux. It is consistent with 
the recent numerical simulations by Komissarov|l7j . 

Let us consider a simple example where the non force- 
free region can be approximated as a thin disk, two di- 
mensional disk on the equatorial plane, dominated by 
the inertia. Then we can suppose no-vanishing B"^ up to 
the equatorial plane with discontinuity at the plane. In 
this boundary value problem, the surface charge, Ue, and 
the current density, , can be deflned on the equatorial 
plane as in the previous section. The energy flux in 9- 
direction is given by Ea. (|33|l . Inside the disk, it is not 
force-free and there is no reason that Bt = vuaB'^ should 
be conserved along the fleld lines and Bt or equivalently 
B'^ develops from zero to flnite value up to the disk sur- 
face. Only beyond the surface where force-free condition 
can be realized, Bt is conserved along the magnetic sur- 
face. The Poynting power measured at inflnity is given 
by 

Pergo = / drrw^FBlB^. (57) 

J Th 

It is interesting to compare the power from the disk 
inside the ergosphere, Pergo, to the Poynting power from 
a black hole Phoie, Ea. (|51|l . As an example we take a 
magnetosphere similar to that suggested by BlandfordQ- 
Using a set of approximations for the magnetic fleld and 
rip for a numerical estimation [s^, it is found that the 
ratio of the Poynting power from the ergosphere to that 
from the black hole, Pergo{o,)/ Phoie{o,), is increasing as 
the angular momentum parameter a increases. The ratio 
of the angular momentum flux, Pergo{o,) / Pl^^^^^d) , is also 
increasing with d. In this simple analysis it is observed 
that for a maximally rotating black hole the power from 
the disk inside the ergosphere can be as much as 30% of 
the power from the black hole. 

Physically Pergo is a part of the gravitational energy 
of the particles of the negative energy orbit in a disk on 
the equatorial plane tapped by the magnetic fleld. For 
example, if the energy momentum tensor of the two di- 
mensional disk is dominated by the inertia, mass density 
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pm, we can make use of Eas. lj^ and We can see 

that the second terms in these equations correspond to 
Pergo and Pergo which depend on the energy(negative) of 
the orbit and the mass accretion rate. The estimation of 
the first terms in Eas. (|38|l and H39|) inside the ergosphere 
is an interesting subject to be studied in the future to get 
a more reahstic estimation of Pergo- 
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